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Describing the Flow of Holes and Electrons in Semiconductors 
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(Manuscript Received June 22, 1951) 

The subject equations are investigated with the aim of establishing some 
general properties of the flow fields which they describe, including the existence 
or non-existence of classes of exact solutions having certain formal properties. 
The results include a number of geometric characteristics of the vector fields 
involved, a suggestive reformulation of the partial differential equations re- 
stricting carrier concentration and electrostatic potential, and several classes of 
exact solutions involving arbitrary constants and/or functions. Of particular 
interest is a family of solutions in closed form for the steady-slate, no-recom- 
bination case involving an arbitrary harmonic function in three dimensions. 
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A. Introduction 

THIS paper is concerned with the system of relations describing the flow 
of holes and electrons in the interior of a homogeneous semiconductor 
subject to the assumption of constant temperature, electrical neutrality, 
and constant difference in concentrations of ionized donor and acceptor 
centers. These relations are: 

div|| P = -eU+^1 (1) 



div IU = e [<R + |2] (2) 

1174 
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° kT 1 

||p - -» p e Ugrad <0 + — grad />J (3) 

° &r 1 

||» = -^ e n grad 1) — — grad n (4) 

» — p = n — Po = N (a constant) (5) 

n, p > (6) 

o o o 

II = Up + II- (7) 

wherein 
n: concentration of negative carriers (electrons) 
p: concentration of positive carriers (holes) 
no', thermal equilibrium value of n 
p : thermal equilibrium value of p 

o 

\\ p : hole current density vector 

o 

||„: electron current density vector 

o 

|| : total current density vector 

/: time variable 

c: magnitude of electronic charge 

k: Boltzmann's constant 

H P : hole mobility constant 

H„: electron mobility constant 

T: absolute temperature (assumed constant with time and uniform) 

V: potential of electrical intensity field 

(R: electron-hole recombination rate function (will usually be regarded 
as depending on p — p and n — n or equivalent variables). 

These relations have fundamental application to transistor electronics, 
photoelectric effects, and related phenomena. Detailed discussions of their 
physical bases will be found in References 1 and 3. In brief, (1) and (2) 
are conservation conditions for the positive and negative carriers; (3) and 
(4) express the dependence of the local current densities on the electrostatic 
potential gradient and on the carrier concentration gradients (i.e., on con- 
duction and diffusion); (5) expresses the condition of electrical neutrality 
under the assumption of a constant difference in concentrations of ionized 
donor and acceptor centers; and (6) and (7) are self evident. 

The present study is directed toward the discovery of (1) general proper- 
tics of the flow fields inside semiconductors and (2) families of exact solu- 
tions to the flow equations. The approach to the latter objective is through 
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the "inverse method" which has proved very useful in the study of vari- 
ous non-linear partial differential equation systems in mechanics. In the 
inverse method, one proceeds by formal devices suggested by the equations 
under study to try to find families of solutions to the equations which in- 
volve arbitrary constants or, preferably, arbitrary functions. This is done 
without reference to any preconceived boundary value problems. After a 
pool of such families of solutions is available, it can be examined from the 
point of view of finding boundary value problems of interest consistent 
with any of the solutions in hand. The likelihood of finding solutions of 
interest in this way is of course greatly enhanced when the solutions in- 
volve arbitrary functions. Aside from providing solutions of some useful 
boundary value problems, the solutions found by the inverse method 
constitute a reference bank of non-trivial exact solutions against which 
to check numerical methods and approximation schemes (based, for ex- 
ample, on the assumption that a particular term can be neglected) for 
solving problems of more immediate practical interest. 

J. Bardeen has demonstrated (in Reference 2) how the steady-state be- 
havior of contact-semiconductor combinations can be explained on the 
basis of the characteristics of (1) the flow field inside the semiconductor 
and (2) those of the barrier layer at the contact. The present study is con- 
cerned in this connection only with the first of these influences. It provides, 
for example, a complete solution for the spherically symmetric flow field 
without recombination for arbitrary currents — a generalization of the zero- 
total current solution given by Bardeen. In the absence of surface recom- 
bination this spherically symmetric solution provides the hemispherically 
symmetric flow field in the neighborhood of a point contact on a plane 
surface and remote from other electrodes or surfaces. This spherically sym- 
metric solution is contained as a particular case in a family of solutions 
involving an arbitrary harmonic function in three dimensions. Other choices 
of the harmonic function can be made to yield flow fields associated with 
numerous electrode configurations of immediate practical interest, for ex- 
ample that of the type-A transistor. 

The objective of the present paper is to find (or establish the non-exist- 
ence of) broad classes of solutions, and not to undertake detailed studies of 
any particular solutions. Such detailed studies of particular cases from the 
family of solutions mentioned above (and from other families found in 
this study) will form the subject matter of papers dealing with specific flow 
field configurations. However, in order to illustrate the interpretation of 
mathematical arbitrary constants in terms of basic physical parameters, 
the analysis of the spherically symmetric solution mentioned above is car- 
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ried up to the point of actual substitution of numerical values in the 

formulae. 

Note: In the following, functions and constants described as "arbitrary" 
are to be considered as being subject nevertheless to the restrictions 
implied by (6). In any particular case it is an elementary matter to 
determine these restrictions and we shall not usually carry out this 
detail. Also, "arbitrary" functions are subject to appropriate dif- 
ferentiability conditions readily evident in any particular case. 

B. Some Properties of the Current Density Vector Fields 

o o 

Several interesting properties of the current density vector fields || p ,|| „ , 

o 

and 1 1 are easily found from (3)-(5). 
It is evident that (3) and (4) can be rewritten as 

\\ p = -en v p grad \V + -j In pj (8) 



and 



| n = — e/j.„ n grad ( V — — In n ) . (9) 



From (3), (4), and (7) we have 

o 

|| = — e(ji n n + n P p) grad V + kT grad (p n n — n P p) (10) 

which because of (5) can be rewritten as 



= — e (ji n n + m p p) grad 



V- — ^_A» in G* » + Mp f) 
e n n + Up 



• (ID 



Now (8), (9) and (11) are all of the form 

U = grad \p 
and hence obviously satisfy the condition 

U • curl u = 0. 
Therefore we have 

CO o 

Theorem 1: \\ p , \\ n , and || are surface-normal vector fields. 
From (8)-(10) we find, using (5) 

o 

curl |j p = — eupgrsLdp X gradl), (12) 
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curl | | w = -en n grad p X gradt), (13) 

o 

curl || = -e(ji n + Up) grad p X grad V, (14) 



and 



whence 
Theorem 2: 



curl ||p curl ||„ curl 



o o o 

|_That is, curl || p , curl || n , and curl || are constant multiples of one another, 
and 

o o o 

Theorem 3: || p , ||„ , and 1 1 are irrotational if and only if 

grad p = (p = p{t)) 

or grad V = (V = V(t)) 

or 13 = 1)(fc /). 

The following interesting relations can be obtained from (8) and (9) 
(they are really consequences of Theorem 1): 

o o 

curl ||p = grad In p X || P (15 

and 

o o 

curl ||„ = grad In n X || n . (16) 
Now from (3) — (5) we find 

|°|p X |°|„ = en n n P kT(n + p) grad p X grad V (17a) 

= y n n P kT(n + p) grad (n + p) X grad V (17b) 

= ien„n P kT grad (» + p) 2 X grad D (17c) 

= \enn».pkT curl [(» + p) 2 grad D] (17d) 



and 



and 



\P I in 



= grad 
fi p e fi n e 



" kT 1 

NV - — (n + p) (18) 



— + — =-(«+*) grad U (19) 
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{Note: As is suggested by (18) and (19), the total carrier concentration 

(? = n-\- p = N -\- 2p = In - N (<P > | N |) 

will frequently appear as the "natural" concentration variable in 
the relations with which we shall be working. Hence, expressions 
involving p, or p and n will often be replaced in the sequel by their 
equivalents in terms of the variable (P. It will be noted that 

grad (P = 2 grad p = 2 grad'ra.] 

Equations (17) and (19) yield at once the following theorems: 

Theorem 4: The vector field 

o o o o o o 

II vll — II V II — II vll 

\\p S\ \\n — || /N J I n I I P -^ II 

is solenoidal. 
Theorem 5: The vector field 

o o 

/ ||p _ ||n \ j s j rrotat j ona i w i tn a potential (— eNV + kT6>). 
Theorem 6: The vector field 



— + — ) is surface-normal (to the surfaces of constant V). 

Mp Mn/ 



\_Theorem 7: \\ p , \\„, || , grad V, and grad p are coplanar vectors. 

o o o 

Theorem (V: The flow lines of any two of the fields || p ,||„ , and || coincide 
if and only if 

grad p = (p = p{l)) 

or grad V = (V = V(t)) 

or V = V(p, /). 

Also, from (17) and (19) we obtain the curious relations: 



° ? , , , i \ 

II kT / \ 

\\p x N_" = _ « g r ad <p X ( ii-' + ^ ) 



^r _ , /IL . I In 

= — -— (P curl I — + — 

2 \Mp Mn 



= - T curl 



(Pi US + li! 

L \Mp Mn/J 



(20a) 
(20b) 
(20c) 
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Finally, by taking the divergence of (7) and making use first of (1) and 
(2) and then of (5), we obtain: 

o 

[Theorem 9: The vector field || is solenoidal. 

C. Formulation of Partial Differential Equation System 
Restricting (P and V 

A very convenient formulation of the partial differential equations re- 
stricting (P and V is suggested by (18) and (19). Taking the divergence of 
these equations and substituting (1) and (2) into the results we obtain: 

div grad (nv - ^ (?) = - <* (<R + i ~ J (21) 

and 

div (<P grad V) = P U + \ ^ j (22) 

wherein for brevity we have set 

J , 1 
a = 1 

Mp Hn 



and 



Mp Mn 



and shall henceforth assume /3 j£ 0, i.e., \i v j^\i n . Equations (21) and (22) 
yield immediately a derived equation not containing explicitly the terms 
introduced by recombination and time variations: 

div grad f NV - — (?) = - ^ div (<P grad V) (23a) 

or 

div \(n + | (Pj grad V - —grad (? \ = (23b) 



or 



div ( ( (P + — ) grad 



v-W^L + w 

oce \ a 



= 0. (23c) 



Either the set (21) and (22) or one of the forms of (23) together with either 
(21) or (22) constitutes a basic set of two partial differential equations deter- 
mining (P and V. We are here considering 61 as (?((?). 
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It will be observed that (23) is equivalent to the condition 

div f| = (24) 

established as Theorem 9. 
(In terms of (P, (10) becomes 

1^ 17* <P + n) grad V - ^ grad (P 



I i - ir t iv i grau u — — gra.u ir i . (25) 

In most of the following sections we shall find it expedient to consider 
separately the cases N ^ and N = (associated respectively with semi- 
conductors of the extrinsic and intrinsic conductivity types). For the case 
N 7^ 0, use will be made frequently of new dependent variables 11 and 30. 
defined by: 

1L^^(P (26) 

eN 



kT 

Kau- TT <P = *0 - 11. (27) 

eN 



That is, 



»-g» (28) 

*U = <U + 3C (29) 

will be substituted into relations involving (P and V to obtain the corre- 
sponding relations in terms of 11 and 3C. Incidentally, it will be noted that 
It and 30 have the dimensions of voltage. 

In terms of 11 and 3C the basic equations (21)— (23) can be written: 

*-« --»[• + »£] (30) 

div hi grad CI + 10] . ?*£ [« + H™] (31) 

div grad 3C + — It grad (it + 0C) = (32) 

wherein (R will be considered as (R(1t). 

It will be observed that, in the absence of recombination and time varia- 
tion, (30)-(32) reduce to 

div grad 30 = (33) 

and [N 9* 0] 

div [It grad (11 + 30)] = (34) 
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The elegant form of this set of equations furnished the original motivation 
for the introduction of the variables 11 and 3C. The comparable equations 
for N = are 

div grad (P = (35) 

[N = 0] 

div [(P grad V] = 0. (36) 

D. The Recombination Rate Function (R 

In order to avoid undue confusion in the sequel we shall at this point 
make some clarifying remarks concerning the function (R. As was stated in 
the Introduction, we basically regard (Rasa function of p — pa and n — n Q . 
However, because of (5), any expression in p — p and n — n can be re- 
placed by one in which (say) p is the only field variable quantity. It is then 
convenient to regard (R as a function of p and write it (R(p). When dealing 
with expressions in terms of (P and of 11, it is convenient to regard 01 as a 
function of one of these variables and to indicate this fact by writing (R((P) 
or 61(11). When we do this we do not mean that (R((P) (say) is the same 
algebraic function of (P as (R(p) is of p, but rather that (R(p) is the function 
of p obtained when one substitutes (P = N + 2p into (R((P). 

For example, for constant mean lifetime recombination 

Oi(p) =-(p- #>) (37a) 

To 

(R((P) ^ i ((P - (Po) (37b) 

<R(11) = ^ (11 - Ho) (37c) 

with to constant; 

and for mass-action recombination 

(R(p) = — [pip + N) - Poth] (38a) 

«0T0 

«<»> = 2r,«Po'+ N) (S? - «* (38W 

m(<a) s 2kVr!Z + N) ^ ~ %t >- (38c) 
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E. Addition of Arbitrary Time Functions to V and 3C 

Since only the gradient of V appears in the basic equations (21) and (22), 
it is evident that if 

"0 = V(x, y, z, t) 

and 

(P = (P(x, y, z, t) 

are a pair of functions satisfying (21) and (22), then so also are 

=0 = V(x, y, z, I) + m(t) 

and 

g> = (P(x, y, z, t) 

where m(f) is an arbitrary time function. 
And since V = 1L + OC, if 

3C = 3C(x, y, z, I) 
and 

K = %(x, y, z, I) 

are a pair of functions satisfying (30)-(32), so also are 

3C = 3C(x, y, z, t) + m{t) 

and 

<U = 0l(*, y, z, /). 

These arbitrary additive functions with zero gradients are physically 
trivial in that they merely reflect the arbitrariness of the reference voltage 
level. They will, however, be retained for the sake of formal completeness 
whenever they appear in the subsequent analyses. 

F. Summary of Solutions for No Recombination or Time Variation 

The next ten sections of this paper (Sections G-Q) contain a sequence 
of detailed analyses in which is determined the existence or non-existence of 
solution fields having certain prescribed formal properties. In most of these 
studies time variability and recombination are admitted and the analysis 
includes the establishment of the class of recombination rate functions 
(R consistent with the property under consideration. In those cases where 
solutions are found to exist, they are expressed in the simplest convenient 
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terms: in closed form, or as solutions of an ordinary differential equation, 
or as solutions of a single partial differential equation. The solutions found 
usually involve arbitrary constants and/or arbitrary functions of various 
kinds. 

The present section is intended to provide a skimpy but compact sam- 
pling of the results obtained in the next ten sections. It will be confined to 
a simple listing of solutions found and furthermore will contain only the 
forms to which these solutions reduce when recombination and time varia- 
tion are excluded. (Some solutions are lost under this reduction.) A heading 
will indicate the section (s) from which the solution comes as well as the 
formal property associated with each solution. 

For the sake of conciseness and simplicity the symbols denoting arbi- 
trary constants and functions in this section are independent of those em- 
ployed in the later sections. They are to be interpreted as follows: 

A, B: arbitrary constants 

h (x, y, z) : any harmonic function 

(or with subscript) 

(13, (P): any given solution field 



[G. grad 13 = 0] 



[H, I. grad (P = 0] 



[J. grad 3C = 0, N ?± 0] 



[K, L. 13 = 13((P), N ?* 0] 
(A ^0) 



13 = A 
|_(P = h(x, y, z) 



[13 = h{x, y, 
|_fl> = A 



z) 



V = A -{- \/h(x, y, z) 



Ne 
(P = ^\/h(x,y,z) 



it \ \ a >.[ B ~ h(x, y, z) 
V = h{x, y, z) + ^A 4 



v = x±AK\ B - h ^y> z) 

kT A 
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(For definition of function A see Equation (87) and Figs. 1 and 2.) 
[K, M. V = V(<P), N = 0] 

Pi) = A In h(x, y,z) + B 

JP = h(x, y, z) 

[N, O. grad (P-grad 1) = 0] 

"V = h{x, y, z) 

(P = h{x, y, z) 

provided 

Lgrad hi(x, y, z)-grad fhix, y, z) = ( 

[N. grad <U • grad OC = 0, N ^ 0] 

~V = Vfh(x, y, z) + h 2 (x, y, z) 

jy e 

® = kf ^ lh{x> % Z) 

provided 

|_grad Jh(x, y, z) -grad h 2 (x, y, z) =0 

[P. grad (P-grad h - 0] 

~V = V + h(x, y, z) 

(P = (P 

provided 

[_grad (P-grad h(x, y, z) = 0. 

G. Solutions with V = V(t) 

Our point of view in general is that (P and V (or ^ and 3C) are functions 
of three space coordinates and time, so that V = V(t) implies for example 

that — = — = — =0. That is to say, we now seek solutions for which 
ax dy dz 

everywhere 

grad 1) = 0. (39) 

From (21) and (22) this condition gives us the following restrictions on 
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(P (and noneonT>(/)): 

div grad (P = (40) 

and 

(R(CP)+i^ = 0. (41) 

I at 

By operating with div grad on (41) we obtain 

(R"((P) = 

(we consistently use primes to denote differentiation with respect to the 
argument of a function of a single variable — e.g., 

d(P- 
whence, 

2<R((P) = A(? + B (42) 

(A, B: arbitrary constants). Substituting (42) into (41) we obtain 

52 + ap--* 

at 

whence 

(P = c (x, y, z)i~ At - B/A (A * 0) (43a) 

or 

<P = c(x, y, z) - Bl (A = 0). (43b) 

From (36) it follows that 

div grad c(x, y, z) = 0, (44) 

that is, c must be harmonic. 

In brief, if (R((P) is of the form given in (42), any V(l) and (43) constitute 
solutions to the flow equations for any harmonic c(x, y, z). Other forms of 
(R((P) admit no solutions with V = V(t). 

It is evident that when recombination is absent time variation is also 
absent, and vice versa. The solutions reduce in this case to: 

V = C (C: arbitrary constant) (45) 

(P = c(x, y, z). (46) 
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H. Solutions with (P = (P(/), N^O 



The condition 



grad(P 


= 


(47) 


yields from (21) and (22) 




(iV + - (?) div grad 13 - 


(48) 


and 




(P div grad V = 


[«» + &] 


(49) 


Two cases arise for (ft ^ 0: 




Case /; 







(p = - 



/3-V 



(50) 



and 



Case £: 



div grad 13 = — — (R((P) 



(P = - 



/fcV 



•W+JS-o 



(51) 



or 



and 



D 



— / 



d(P 

2(R((P) 



(Z): arbitrary constant) 



(52) 



(53) 



div grad V = 0. 

When recombination is absent, these cases reduce to: 

(P = £ (£: arbitrary constant) (54) 

and (53). 

When time variation is absent, Case 2 again yields (53) and (54). 

Tt should be recalled that V can depend on / as well as x, y, z; so that 
arbitrary functions of / play the role of arbitrary constants in (51) and (53), 
whenever time variation is allowed. 
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I. Solutions with P = P(t), N = 

For N = 0, only Case 2 of the previous section occurs, because the con- 
dition (P = (implying no carriers!) is of no interest. 

J. Solutions with 3C = 3C(/), N 9^ 
For grad 3C = 0, (30) and (32) yield: 

.» + £=- 

and 

div grad at 2 = 2 div 11 grad 01 = 0. (56) 

Taking the div grad of (55) multiplied by 01 we obtain 

div grad 0101(01) = 
whence, because of (56) 

— 0l(R(0l) = F01 2 + G (F,G: arbitrary constants) 

eN 

or 

^ (R(0l) - Ml + GOt" 1 . (57) 

eN 

Substituting this permitted form for the recombination rate function into 
(55) we obtain 

^ + F01 2 = -G (58) 

at 

whence 

01 = Vf(x, y, z)t~ n - G/F (F 9* 0) (59a) 

or 

01 = V/(x, y, z) - Gt (F = 0). (59b) 

From (56), f(x, y, z) is subject to 

div grad fix, y, z) =0. (60) 

In summary, if and only if Gt(0l) has the form (57), there are solutions 
for which 3C = 3C (t) (arbitrary). The 01 is given by (59) in which/ is an 
arbitrary harmonic function of x, y, z. 
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In terms of (P and V these solutions are given by: 

eN 
(P = j~ V/(x, y, «)«""- G/F (F * 0) (62a) 



or 



and 



or 



eN 



9 = kf Vf{Xi y ' Z) ~ Gt {F = 0) ' (62b) 



V = 3C(/) + y/f{x, y, z) € -" - G/F (F ^ 0) (63a) 



V = K(t) + Vf(x, y, z) -Gt (F = 0). (63b) 

For no recombination ((R = 0), these results specialize to (59b), (60), 
(62b), and (63b) with G set equal to zero. It should be noted (see (55)) that 
absence of time variation implies absence also of recombination. 

K. Solutions withD = V((P, t), grad <P ^ 

In Theorems 3 and 8 of Section B we have shown that some very inter- 
esting properties are implied by the condition 

grad V X grad (P = 0. (64) 

In sections G-I we have treated the cases grad V = and grad (P = 0. 
We now turn to the remaining possibility leading to (64): 

V = V(<P, t) with grad (P ^ 0. (65) 

Substitution of (65) into (23b) leads to 

Two cases arise. 
Case 1: 

jr + S^g-S-a 

8 /d(P e 
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This condition clearly satisfies (66) and leads to 

(g(/): arbitrary function). (67) 



1)((P, /) = g(0 + ^ In 
ae 



<p + ** 



The restriction on (P is then provided by the result of substituting (67) 
into (21): 



[• 



div grad I (P — — In 
a 



(P + — 



H^+ig. 



(68) 



Any (P satisfying (68) constitutes with (67) a solution having the property 
desired. 

If (65) is substituted into (25) it will be found that the condition 

T 8 / d<P e 

o 

is equivalent to || = 0, so that Case 1 is characterized by zero total current. 
Case 2: 



|3 /3(P e 



0. 



In this case (66) can be written in the form 



div grad (P _ d . 

(grad (P) 2 " d(? 



8 J d(P e 



= +(<?, I). (69) 



From (69) it follows that (P must be of the form (P (/?, /) with 

div grad h(x, y, z, t) = 0. (70) 

In summary we have 

o 

V Theorem 10: If V = T)((P, /) with grad <P ^ 0, then either || = or V = 
|_ V(h, I) and (P = (P(/i, /) with div grad h{x, y, z, I) = 0. 

We shall investigate the restrictions on the functions h(x, y, z, /), V(h, /), 
and (?{h, l) in the next two sections. 

Theorem 10 remains unchanged if recombination is absent. If time varia- 
tion is absent, it simply drops / as a variable in the functions mentioned in 
the theorem. If both recombination and time variation are absent, the 
theorem can be strengthened to: 

Theorem 11: If both recombination and time variation are absent and 
V = V(G>), then V = V(h) and (P = <?(h) with div grad h(x, y, z) = 0. 



[ 
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L. Solutions with V = V(h, I), (P = (P(h, I), grad (P j* 0, 

DIV GRAD h = 0, N 9^ 

For formal reasons we shall work, not with the conditions (P = (P(/z, t) 
and V = V(h, /), but with the equivalent conditions 

01 = Ol(/z, /) and 3C = 3C(A, /). 



(71) 



The condition grad (P ^ now implies -=- ?^ 0. 

o/i 

Substitution of (79) into (30) and (32) yields— after use of (70): 



— (grad ft) - - - 



0t(ll) 



2*r a/ 



eiv aota/* 
2*raA a/ 



and 



a 
dh 



' pkT 
ae 



+ 01 



5+»S>- 



(72) 



(73) 



From (73) we get 



d_3C 
dh 



**-** 



fikT 



+ 01 



(74) 
(j(t) : arbitrary function) 



which yields upon substitution into (72) : 
(grad hf 



d 

dh 



*«-«S 






+ 01 



(75) 



in which 01, — r , and — — are, of course, functions of h and of /. 
dh dh- 

in determining the combined implications of (75) and (70) three cases arise 

r)-TC 

according to whether or not -—-^ = or grad (grad h) 2 = 0. 

on 

Case 1: 



a// 2 * ' 



grad (grad /z) ^ 0. 
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In this case no satisfactory interpretation has been found when time 
variability is present. 
When time variation is absent, we work with the conditions 

01 = 11(A); JC - 3C(h) 

with 

div grad h{x, y, z) = (76) 

and arrive at counterparts of (74) and (75) : 

3C' . H ~ 3J3L ( y = Pkf) (H: arbitrary constant) (77) 

7 + %■ \ « e / 

and 

3C" (grad hf = ( g 7 ^' )' (grad hf = - « <R(«U). (78) 

From (78) it is evident that (R ^ implies 3C" ?^ and grad h ^ 0. So 
we have 



(grad A)' = /0 ff ., Mtx . (79) 



/ S - w \ 
W + ii ) 

which is of the form 

[grad A(s, y, z)] 2 = 4>{h). (79a) 

Now from (79a) follows 

grad h X grad (grad hf = (80) 

which implies that the vector lines of the field grad h are all straight. Since 
h is harmonic, this restricts the choice of h to the potential fields associated 
with a uniform parallel flow, a straight line source, or a point source. Hence, 
for suitably chosen rectangular coordinates (x, y, z), circular cylindrical co- 
ordinates (p, 0, z) or spherical polar coordinates (r, 6, <f>), the only possibili- 
ties, are respectively 

h = x -> (grad hY = 1 (81a) 



or 



h = In i -* (grad A) 2 = -, = ™ (81b) 

P P 
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or 

h = - -> (grad hf = - 4 = h\ (81c) 

The possibility h = x violates one defining condition for the present case 
(i.e., grad (grad h)- ^ 0) and hence will be left for consideration in Case 3. 
The remaining two possibilities lead respectively to the following forms of 
ordinary differential equation for the determination of 11(A): 

(Trr)' + r-««- (82b) 

or 

Given any H(/j) satisfying one of these equations, the associated 3C(h) is 
obtained by integration from (77): 

3C(A) = f( H ~/^f ) dk + J (J ... . .. (83) 

J \ 7 + 1L / (/: arbitrary constant). 

It is evident from (72) that Case 1 does not exist if both recombination 
and time variation are absent. 
Case 2: 

dh 2 

r)TC 
In considering this case we shall exclude the condition — = because it 

on 

has been included in Section J. 

From the condition — — = we have 
air 

3C = k(t)h + C(l) (k(t), /(/): arbitrary functions) (84) 

with k ?± 0. This shows that 3C itself is a harmonic function and we can with- 
out loss of generality use it in place of h. 

Equations (74) and (75) now yield the two conditions on 11(30, /), (R(H), 
and.TC(.r, y, z, /): 

■/.\ n dIL 

>®-*ax = . (85) 

11 + 7 
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and 



<*<*) + ^ (***!? + ^) = 0. (86) 



4 

3 
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Fig. 1 — The transcendental function A (x). 
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Fig. 2 — The transcendental function A (.r). 

For the integration of (85) we need the transcendental algebraic function 
of a single real variable defined by 

A(.r) + In | A(x) - 1 | = x. (87) 

This function is plotted in Figs. 1 and 2. It will be observed that x is always 
a single-valued function of A; while A is a single- valued function of x for 
x > 0, a double-valued function for x = 0, and a triple-valued function for 
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x < 0. The single-valued monotone functions Ai , A 2 , and A 3 are denned 
respectively by the restrictions A > 1, 1 > A > 0, and A < 0. When A is 
used without subscript it is implied that either A x , A 2 , or A 3 can be used. 
It will be useful to remember that 

A'W = A -^r^. (88) 

A(.v) 



In terms of the function A, (85) integrates to 



01(30,/) = [j(t) - y] A 



m 



(/) - 3C" 



(89a) 



_ j(0 - 7 _ 0" 5* 7) 

(m(t) : arbitrary function) 

or 

<U(3C, = m{t) -X. (i = 7). (89b) 

The latter case (j - y) corresponds to V = V(t) and hence was included 
in Section G. Therefore in the following we shall consider on\y j ^ y. 
Now by making use of (89a) and (88), (86) can be rewritten in the form: 

2kT 0161(01) . j'OL 2 d3C,/(w-3C)_ , 

_ — . -+- — — ' — r~ ■ — 7/1 

eN 01 -j + 7 0'- 7)CU- -j + y) dt 7-t (gQ) 



f primes denoting here — 1 



We now observe that the right side of (90) is harmonic, while the left side 
is a function only of 3C and /. From this it follows that the right side can be 
written in the form: 



ax /'( fit - 3C) _ , _ _, A [tn - oe" 

~T~~ — ?rl 



' " < (,) [j 



dt j - y L J - 7 _ 

From (90), (91) and (89a) follows 

^(R(«U) - —J-% 

eN j — y 



+ r(t). (91) 



(92) 



+ 0i-i + 7/ +T J_ ni + glnpL.-! 

01 \ ] - y \j - y 

Since (92) is of the form 

(R(0l) = 0(11, 0, 

The result of taking ( — J of the right side must be zero identi- 



It = const 
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cally in 11. Making use of the algebraic lemma that 



A + Bx + - x + (d+- x ) In 



»- 



= 



implies A=B = C = D = E = 0, we arrive at two possibilities: 
Possibility 1: 

(q(l) = r(t) = and j - y = K<T U ) 

(K, L: arbitrary constants) 
This yields from (92), (89a) and (91): 



^ 01(01) = ill 

eN 



11(00, 1) = Ki 



■Ai 



(m(l) - 3C) 



] 



(93) 
(94) 



and 



3C(*, y, z, t) = f" s(x, y, z) + C hi f%" Wrntt) + m'{t) dt (95) 

where s(x, y, z) is any harmonic function. 
Possibility, 2: 

(j(t) = M,q(t) = Q,r(t) = R) 

(M, Q, R: arbitrary constants) 
This yields from (92), (89a) and (91): 



eN 1L 



R-\ 
11(30, /) = {M - 7) A 



M — 7 


+ <21n 


~m(t) - 30.' 
_ M — 7 





It 



M — 7 



and 



K(x,y,z,t) =e QmM - y) u(x,y,z) 



■ Qtl(if-y) I -(.Qtl(M-y)) 



R + w'(/) + 



Q 



M — 7 



;« 



(0 



<//. 



(96) 



(97) 



(98) 



where u(x, y, z) is any harmonic function. 

In the absence of recombination, Possibilities 1 and 2 lead to the same 
result: Equation (97) and 



JC(.t, y, z, t) = u(x, y, 2) + m(l). 



(99) 
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In the absence of time variation, (86) shows that recombination is neces- 
sarily absent, too, so the results reduce to 

1l(3C) = AA (?-Z^) (100) 



with 3C(x, y, z) any harmonic function and A and B arbitrary constants. 
This solution for the case grad 50. 5* 0, together with that given by (59b) 
and (60) (with G = 0) for the case grad 3C = 0, constitute a veritable gold 
mine of useful solutions because of the arbitrary harmonic function involved. 
An example involving a particular choice of 3C will be examined in Section R. 



Case 3: 



,2, 



|p * 0, grad (grad h) 2 = 0. 



In this case (grad h) 2 is a function of / so that (75) can be written in the 
form 

From this it follows (because div grad h = 0) that 

h{x, y, z, I) = d{t)b(x, y, z) + l(t) (101) 

with 

div grad b(x, y, z) = 0. (102) 

The condition grad (grad h)"- = now requires further that 

grad (grad by = 0. (103) 

But any 5(.v, y, z) satisfying both (100) and (101) can, by suitable choice 
of axes, be written 

b = Sx (S: constant). 

This leaves us with exactly the same totality of solutions as we could have 
obtained by setting It = 11 (.v, /), 3C = 3C(x, i) in the first place. So we replace 
h by x in (74) and (75) and obtain: 

•M - — 
33C = m _ dx (104) 

dx 7 + 11 
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and 



dx 



dx 

L 7 + 11 J 






(105) 



Any llfo /) satisfying (105) can be substituted into (104) to obtain 3C(x, t) 
from 



3C(*, t) = /(/) + f 






</a: 



(106) 



7 4- <U 

(/(*): arbitrary function). 



If recombination is absent, (R(1l) disappears from (105). If time variation 

is absent, — disappears from (103) and j(t) and/(/) are replaced by arbitrary 
dt 

constants. In the latter case, the standard change of variables 



WCll) for 



dx 



i \ d d 

^M for Tx 



(107) 



reduces the solution of the second order equation (105) to the solution of a 
first-order equation followed by a quadrature. If both recombination and 
time variation are absent, the substitution (107) reduces the solution of (105) 
to two quadratures. 

A set of equations equivalent to the steady-state ( — = J forms of (104) 

and (105) has been the subject of an extensive numerical investigation by 
W. van Roosbroeck (Reference 1) for the recombination rate functions 
given in (37) and (38). 

M. Solutions with V = V(h, /), (P = (P(h, /), Orad (P ^ 0, 
Div Grad h = 0, N = 



For these conditions (21) and (23b) yield 



a 2 (P f . M 2 oce 

W- (grad h) = kf 



,1 fd(P dh , d(P\~] 



2 \dh dt 



dt 



and 



d_ |~d<P _ <xe dV 
dh \_dh $kT dh 



(grad 



hf = 0. 



(107) 



(108) 
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Since we do not here allow grad h = 0, (108) implies 

dV y [dh~ |W J PkT fm/A ... . .. v (109) 

— = — t =■ 7 = - — (|(0 : arbitrary function) 

d/j (P ae 



Case i: 



2 

^ ^ 0, grad (grad hf * 0. 



In this case, as in the associated case in Section L, the implications of 
(107) together with 

div grad h(x, y, z, t) = 

are not known when time variation is present. 

When time variation is absent, we work with the conditions 

(P = (P(A) and V = V(h) 

with 

div grad h(x, y, z) = 

and arrive at counterparts of (107) and (108): 

fl>". (grad hf = %(R«P) (110) 



1 -_./\/ 



and 

(<$>' -i(TO'J' = 0. (Ill) 

Proceeding as in the analysis of Case 1 of Section L, we infer that // must 
be of the kind given by (81b) or (81c). The associated second-order differ- 
ential equations restricting <y(h) are then, respectively: 



and 



<P"-^«-*<R((P) =0 (112a) 

kT 



(P"-~(H((P) =0. (112b) 

kT h* 



The V(h) associated with any solution of (112) can be obtained by integra- 
tion from 

V(h) = C + [ y6> ~ — dh (C, D: arbitrary constants). (113) 



1200 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1951 

It will be noted from (110) that simultaneous absence of recombination 
and time variation is inconsistent with the defining conditions of this case. 

Case 2: 

dh 2 

We shall exclude the possibility of — - = because it is included in Section 

oh 

I. Then, proceeding as in Case 2 of Section L, we conclude that (P itself is 

a harmonic function and can be used in place of //. (107) and (109) 

then become- 



and 

dV y[l - |(01 



(115) 



d(p (p 

or 

1)((P, /) - t[1 - 1(0] hi (P + }(0 0(0: arbitrary function). (116) 

d(P 
Because — is harmonic and a function of (P, we have 
dt 

2(R((P) = - d -^. = E(P - P 
dt 

(E, F: arbitrary constants) 
2(R((P) = E(P - P (117) 



whence 



and 



or 



(P(x, "v, z, I) = e-*'m(x, y, s) - L (£ ^ 0) (118a) 

E 



<?(x, y, 2, = m(x, y, z) + Ft (E = 0) (118b) 

where tn(x, y, z) is an arbitrary harmonic function. 

If recombination is absent, these results specialize to (116) and (118b) 
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with F — 0. If time variation is absent it follows from (114) that recom- 
bination is absent, too, and the results specialize to 



1)(<P) = G + 6 In <P 



(119) 



with (P(x, y, z) any harmonic function and G and R arbitrary constants. 
These solutions play the same role for the intrinsic semiconductor (N = 0) 
that (100) does for the extrinsic (N ^ 0). 
Case 3: 



d 2 (P 
dh 2 



^ 0, grad (grad h) 2 = 0. 



In this case it can be shown, just as in Case 3 of Section L, that no gen- 
erality is lost by considering (P = (9{x, t) and V = V(x, t) in place of (P(h, t) 
and V{h, /). Equations (107) and (109) then become 



and 



d 2 (P _ ae 
dx* ~ kf 


[«« + &] 


dv _ y 


K-H 




dx 


(? 





(120) 



(121) 



Any solution of (120) when substituted into (121) gives an associated V 
from 



d(P -, A 
V(x, t) = q(t) + 7 / ^ 



rfa;. 



(122) 



If recombination is absent, (R((P) merely vanishes from (120). If time varia- 
tion is absent, the functions g(/) and q(/) are replaced by arbitrary constants 
and the standard change of variables 



u((P) for 



d(P 
dx~ 



U((P) ^ for £ 



(123) 



leads to a solution of (120) in two quadratures. An equivalent solution is 
given by W. van Roosbroeck in Reference 1. From (120) it follows that 
recombination and time variation cannot simultaneously be absent for 
Case 3. 



1202 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1951 

N. Construction of Solutions from Orthogonal Harmonic Fields, 

N * 

There are many known examples of pairs of harmonic functions h x (x, y, z) 
and h 2 (x, y, z) that have orthogonal vector fields— that is, for which 

grad Ai ■ grad fe = (124) 

with grad h ^ and grad h 2 9* 0. [E.g., the real and imaginary parts of any 
analytic function of a complex variable.] From any such pair of functions 
we can construct the following solutions of (33) and (34) : 

11 = h ; 3C = fa- h (125) 

and 

11 = Vh ; 3C = h . (126) 

In terms of (P and V these solutions are 

(P = ^h; V = h 2 (127) 

and 

(P = ^VAi; v = Vh + ?h. (128) 

kl 

The validity of the solution (125) is seen from (33) and this expanded 
form of (34): 

% div grad 11 + grad It -grad (it + 3C) = 0. (129) 

Similarly, the validity of (126) follows from (33) together with a different 
expansion of (34): 

div grad It 2 + 2 grad It -grad JC = 0. (130) 

It is evident that a given h x and h 2 can be interchanged in the above 
solutions to yield different solutions, and also that any given hi or h 2 can be 
replaced by an arbitrary constant multiple of itself plus a second arbitrary 
constant. 

O. Construction of Solutions from Orthogonal Harmonic Fields, 

N = 

We can write the differential equation system for the intrinsic semi- 
conductor [(35) and (36)] in the form: 

div grad (P = (131) 

(P div grad V + grad (P-grad "0 = 0. (132) 
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From these we verify the solution: 

(P = hi ; V = h (133) 

for any harmonic hi and h 2 satisfying (124). 
The solutions given by (127) and (133) have the property 

grad (9-gra.dV = 

and so may be considered, in a sense, complementary to the solutions in 
Sections L and M for which 

grad (P X grad V = 0. 

P. Superposition of a Harmonic X Field, N. j± 

Inspection of the equation system [(33), (130)] reveals the following 
superposition theorem for obtaining new solutions from some known solu- 
tions/or the case of no recombination or lime variation: 

[Theorem 12: If [ C U, 3C] is a known solution and if h is any harmonic func- 
tion such that grad 'U-grad h = 0, then [11, 5C + h] is also a solution. 

Or, in terms of (P and V : 

[Theorem 12': If [(P, V] is a known solution and if h is any harmonic function 
such that grad (P- grad h = 0, then [(P, V + h] is also a solution. 

In the latter form it is evident from Section O that the theorem holds also 
for N = 0, but does not extend the results of Section O. 

Q. A Partlal Differential Equation in Terms of 3C Alone, N ^ 

For N = 0, (21) provides a differential equation involving only one de- 
pendent variable — (P. We shall now derive an analogous — but vastly more 

complicated — differential equation for the case N 5^ 0, — = 0. 

dt 

For this case (30) and (32) become 



div grad 3C = — — 01(11) 



and 



div grad 3C + - 11 grad (01 + 3C) = 0, 
or in terms of a familiar vector symbolism 
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V 2 3C = - £ (R(«a) (134) 

and 

V- |V3C + - 1tV(ll + 3C) = 0. (135) 

Now let §(11) be the inverse function to 01(11), i.e. the function such that 

s(<R(m)) ■ it. 

Then from (134) we have 

U=sl'--V a acV (136) 



a / 
Substitution of (136) into (135) yields after some computation 

SSV(V 2 3C) - - (SS" + S' 2 )(VV 2 3C) 2 

(137) 



+ S'V3C-VV 2 0C - ^ (8 + 7)V 2 3C = 



where s'(^) = -77 S(^), etc. 
chf/ 



( N , \ 
S, S', S" are considered as given functions of ( — - V 2 3C 1 . 

The simplest meaningful choice of S is 

•(-^*)-V-(-** B ) +J! < 138 > 

(/, A' : prescribed constants) 

corresponding to constant mean lifetime recombination. For this §, (137) 
specializes to 

j(K - 7v 2 3c)v 2 (v 2 3c) - / 2 (vv 2 oe) 2 

+ /V3C-VV 2 3C - (7 + K - 7V 2 3C)V 2 3C = 0. 

If any 0C can be found satisfying (139), the associated It is given (from 
(136)) by 

11 = 7V 2 3C + A\ 
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R. Sample Application of the Results of Section L: Spherical 
Symmetry, N t^ 

As an example of the solutions included in the results of Section L we 
consider the case of a spherically symmetric field about a point (or spherical) 
source of current. 

We take, as the most general harmonic function having spherical 
symmetry, 

5C = L - + M (140) 

r 

(L, fit: arbitrary constants). 

For the time being we shall assume L ^ and M ^ 0. Then from (100) 
and (28) and (29) we have 

V = M( S - it Z L/r ) + M+± (14.) 



and 



-%*<!±*f&Y 



In terms of V and (P, (3) and (4) can be written 

° -y. v e |~, ... , „ . kT 



and 



((P - N) grad "0 + — grad (P (143) 

- (<P + N) grad V - — grad <P (144) 



which yield upon substitution of (141) and (142): 



j p = l^eNli^A - lHr, (145) 



and 



l-^ z (fr I + 1 )^ 



(146) 



where T\ is the unit radial vector. The total current density is obtained by 
adding (151) and (152): 



= ±eNL 



U. + n P ) ^ A + ( Mn - Mp)] j, rr (147) 
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The currents flowing are obtained from the current densities from the 
relation 

I = ftr 2 || -ri 

where ft is the solid angle (with respect to the origin) within which the 
flow field lies. (If the current source is surrounded by the homogeneous 
semiconductor, ft = 4t; if it lies on a flat surface of a large slab, ft = 2t, 
etc.) So we have 

I p -±QH,eNZ(jLl - ij (148) 

/„ m^Q^drzfl-Z + l) (149) 



(150) 



We shall now obtain expressions for the mathematical parameters B, A, 
L, and M in terms of meaningful physical quantities: I p , I n ,"0« and (?«,. 
(Subscript °° refers to values of variables as r becomes very large.) We 
shall take our reference voltage as the voltage "at infinity" so that V M = 0. 
Setting \/r = in (141) and (142) we obtain 



o-n(i^ + 



M 
\ A / 

and 



«.-j£n 



from which follows (for A 9^ 0) 

B=iln|£ + 

and 



(151) 



eN 
Fro:-. (148) and (149) we readily find (for 1^0): 



M = -~<P«. (152) 



A = kT ^I n + ^ p (153) 

e MpA> — MrJp 
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and 



L = 



fipln — litJ t 



Finally we substitute (152) and (153) into (154) to get 



B = 



kT H P I„ + H n Ip . 
6 Mp-^n — Vnlp 



(N - (PjjX P /n + (N + <?»)pnl t 
N(jl p I n ~ Unlp) 



(154) 



(155) 



Equations (152)— (155) give the desired expressions for M, A, L and B 
in terms of I p ,I n , and (P w for V M = if A ^ and L^O. 
For A = we can repeat the above steps using 



and 



V = B 



<5> = ^(B-M- L/r) 



in place of (141) and (142). The result for L ^ and V* = is 
I p = -\&ii p eNL 
In = |flju„ eNL 
I = %QeN(ji n - n P )L 



with 



B = 
M = - — (P. 



and 



Z, = 



- 2L 



2/„ 



£2/t p eA 7 ft/i n eA 7 



MpZ n — MnZ ? 

O/jp /i n <?A 7 



(156) 
(157) 

(158) 
(159) 
(160) 

(161) 
(162) 

(163) 



It is evident that A = implies V = constant and Hpl„ -\- n,J p = 0. 

The condition L — makes 3C = constant, so we use (62b) and (63b) 
and set 



(164) 



and 



Ne 



(P = ^VR + S/r 



V = Q + VR+ S/r 

(Q, R, S: arbitrary constants). (165) 
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From (143) and (144) we obtain 

/„ = - a -g*S (166) 

I. = - ^ S (167) 



and 



1 = " m (Mn + Mp) ~ s - (168) 



From (164) and (165) we readily obtain for V„ = 0: 



i? = (^(Po Y (169) 



and 



= -^(Poo. (170) 

Ne 



It is evident that Z = corresponds to the case 0C = constant and implies 
Hpl n — nJp = 0. 

The foregoing now provides a formal solution with !)„ = for every 
assignment of values to (P« , ip , and /„ . There remains the question of the 
requirements imposed by the condition 

n,p>0 

which is equivalent to 

(P > | N | . (171) 

This implies first of all that (P«, must be chosen > | N | . 

It is instructive to look first at the case L = 0. Equation (165) shows 
immediately that (171) requires the choice of the positive sign for the radical 
for N > and the negative f or N < to avoid <? K < \ N \ . We further find 
by substitution of (166) and (169) into (165) that (171) requires 



-Ufe7>p(<P»-|iVr)_ 



NI V . (172) 



The bracketed factor is positive. Since we are interested only in non-negative 
values of r, (172) imposes no restriction if I p is zero or not of the same sign 
as N. However, for N and I p of the same sign, (172) establishes an inner 
radius inside which the solution does not satisfy (171). This may be regarded 
as establishing the minimum radius for an inner spherical electrode for 
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prescribed I v and (?„ , or alternatively as limiting the possible choices of 
I p and (Poo for prescribed inner electrode radius. Had we chosen the con- 
stants Q, K and S so as to obtain prescribed values of (P and V at a pre- 
selected electrode radius r , restrictions analogous to (172) on the maximum 
radius would appear. 

For the case A = the restriction analogous to (172) is 



r > - 



I P . (173) 



JlupkTif?*, - \N\) m 

Since the bracketed factor is positive, (173) provides no restriction for 
I p > 0, but for I p < establishes a minimum radius of the kind just dis- 
cussed. 
For L,A 9* 0, the analog of (172) and (173) is 

L/A 

_, (kTV„\ ~ (kT\N\\ (174) 



r > 



A-'(fr-)-A-(^) 

\ NeA / \ NeA / 



where A and L are given by (153) and (154) and A -1 denotes the inverse 
function of A — i.e., 



or 



A-'[A(x)j = x 



A~'(A) = A 4- In | A - 1 



Equation (174) is a minimum radius restriction of the same kind as those 
obtaining for A = 0, and L = 0, but the relationship between the minimum 
radius r and (P,, , I P and /„ is considerably more complicated than in the 
more degenerate cases. 

It will be noted that the relation 

kTV* = K (B - M 



eN A \ A 

(with A, B, M given in terms of (P« , I p , I n by (152), (153), and (154)) deter- 
mines which function (Ai , A2 , or A 3 ) is to be used for A in any given case, 
because any assigned value (^ 0) is taken on by one and only one of 
(Ai , A-. , A 3 ). 

If surface recombination is negligible as well as interior recombination, 
this spherically symmetric solution is of use in the study of "point" con- 
tacts on a plane surface of a semiconductor. [Fig. 3 and Ref. 2.] 

The results of this section can easily be duplicated for any other choice 
of the harmonic function UC to obtain a great variety of specimen solutions. 
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Solutions based on OC's having a single source singularity (such as the ex- 
ample above) will contain four mathematical parameters, and hence will 
permit arbitrary selection (subject to (6)) of the physical parameters, I p , 
I n , (Poo , and V w . However, solutions based on 3C's having more than one 
source singularity will provide only a subset of the possible assignments of 
the physical parameters. For example, the harmonic function associated 
with the electrostatic field produced by two separate point charges each equi- 
distant from two parallel infinite plane conductors provides solutions of 



! / 

i / 







Fig. 3— Point source flow field, useful in connection with point contact theory 



w//////sy///////////w^^ 




Fig. 4 — Two-source flow field between conducting planes, useful in connection with 
Type A transistor theory. 

interest in connection with the type A transistor configuration (Fig. 4). 
However, the family of solutions obtained contains only a five-parameter 
subset of the six-parameter family obtainable by arbitrary assignment of 
I p i , I pi , 7„i , /»s , (Poo , and V M . 

S. Sample Application of the Results of Section M: Spherical 
Symmetry, N = 

We now round out the considerations of Section R by exhibiting the 
related solutions for N = (i.e., the intrinsic semiconductor). 
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In accordance with the results of Section M, we choose for (P the most 
general harmonic function with spherical symmetry: 

(P = A - + B (A, B: arbitrary constants). (175) 



From (119) then, for A ^ 



t 1 



V =* Hlh [Aj+B) I -a 



and from (175), (176), (143), and (144) 



Z,-iQMtlf£+Y 



In = ^Hnel(S-^\ 



I = ^QeA 



- kT 

G"n + f*p)H — (n n — H p ) 



From (177) and (178) we obtain 



A m 



fl n Ip — H P In 



and 



H = 



kT H n I p -\- Pplj. 



e H„I p — Up I n ' 

and from (175) and (176) fori)*, - 0: 

B = (Poo 

and 

G= -I~I\nB = - k _I^I P + ^I n 



In (P. 



(176) 

(177) 
(178) 
(179) 

(180) 

(181) 
(182) 

(183) 



e n n Ip — n P I n 
The condition (P M > | N | = introduces the restriction (for 4^0): 

1 



r > 



.fyipMn^rov 



(Unlp — Up In) 



(184) 



Evidently this implies no real restriction for n„I p — Hpl„ < (i.e., A < 0), 
but introduces a minimum radius — of the same kind we have already dis- 
cussed — when n n Ip — Hpl„ > (i.e., A > 0). 
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For A = 0, (P is constant and, by Section I, V is harmonic. So we set 

(P = (P„ > (185) 

and 

V = C - + D (186) 

r 

and obtain from (143) and (144) 

I p = ifi/ip eC^» (187) 

and 

/„ = $anneC<?« . (188) 

From (187) and (188): 

A = lip = 2/» = U nlp + Mp/n ( lg9 ) 

and from (186) for V m = 0, 

D = 0. (190) 

Evidently vt = is associated with the condition 

Pnlp ~ V-pIn — 0. 

T. Summary List of Symbols 

Coordinate Systems: 

(x, y, z) : ordinary rectangular cartesian coordinates, 
(p, 6, z): ordinary circular cylindrical coordinates. 
(r, 6, </>): ordinary spherical polar coordinates. 
Fi : unit radial vector in (r, 0, 4>). 

t : time variable. 

Physical Variables: 

n : concentration of negative carriers (electrons). 
p : concentration of positive carriers (holes). 
(P: total carrier concentration = n + p. 

<U : s*?(P (N9*0). 

eN 

(R: recombination rate function. 
V: electrostatic potential. 
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kT 

3C: m V - 'U = V - — . <P (tf 3* 0). 

«A 
|| : total current density vector. 

o 

||„: electron current density vector. 

o 

||„: hole current density vector. 

subscript "0": designates thermal equilibrium values. 

subscript "co": designates values "at infinity". 

Physical Constants: 

T: absolute temperature. 

e : magnitude of electronic charge. 

k : Boltzmann's constant. 

/*„: electron mobility constant. 

fj.pi hole mobility constant. 

a : ■ l/n P + l/n„ . 

(3 : ■ 1/m, - 1/Mn • (Assumed ^ 0) 

• - gjl 

A 7 ": = n — pa . 

Other Constants and Functions: 

A,B, ■■■ ,Z ((except /, N, and T)), 
A,B,--,Z, 

A, B, • • • , Z: arbitrary constants 

a,b, ■•• ,z ((except e, h, k, n, p, r, t, x, y,z)), 

S,b, ■ ■ ■ , z: arbitrary functions of variables designated (e.g., j(l)). 

h, hi , h 2 : harmonic functions of variables designated at place of usage. 

A: A(x) is defined by the relation A(.v) -f- In | A(.v) — 1 | = x. 

(See Figs. 1 and 2.) 
S: S(1l) is defined by S[(RCU)] = 01. 
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